Retrieval of multiple spin waves from a weakly excited, metastable atomic ensemble 
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The emission of light from a multiply excited atomic ensemble is examined and it is shown 
how symmetric (spin-wave) and non-symmetric states of excitation radiate into spatially separate 
field modes. This observation has potential application to single photon generation and spin wave 
entanglement, since in the presence of atomic interactions it can result in isolated single photon 
emission into a phase-matched field mode. 
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I. INTRODUCTION 

Single excitation states of matter and light are an im- 
portant resource for quantum communication and com- 
puting protocols. In this paper we consider cold, optically 
thick atomic ensembles as the medium for the creation 
of matter excitations and the source of single photons 
emitted into a well-defined mode [l|. The DLCZ quan- 
tum repeater protocol Q employs atomic ensembles as 
quantum memory elements, interfaced with flying light 
qubits: the detection of a Raman scattered photon her- 
alds the creation of a matter excitation in the form of 
a spin wave, and its coherent storage provides a way to 
realize deterministic single photons through a read-out 
protocol involving feedback Q. To avoid noise due to 
multiple excitations of the medium, weak laser excitation 
is employed. This bottleneck limits the speed of single 
photon production providing motivation for alternative 
approaches, for example, harnessing atomic interactions 
in order to limit the deleterious effects of laser induced 
multiple excitations. 

The Rydberg blockade mechanism, in which the exci- 
tation of a particular atom prevents excitation of nearby 
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atoms in a small sample, is based on strong di 
interactions between atomic Rydberg levels 
subsequent decay of the excited atom will generate a sin- 
gle photon, without the need to rely on repeated weak 
laser excitation cycles. An alternative proposal, operat- 
ing with a larger sample, outside of the Rydberg block- 
ade regime, involves using Rydberg atom resonant dipolc- 
dipole interactions to decouple the decay of multiply ex- 
cited atomic configurations from the preferred output ra- 
diation mode, leaving only the singly excited component 
to decay with the emission of a single photon [9(. In 
this paper we discuss this radiative decay mechanism in 
detail, showing explicitly how interaction induced atom- 
pair phase shifts decouple multiply excited atomic states 
from the phase-matched mode. As a consequence, we 
show that the quantum statistics of the phase matched 
field mode has a strongly single-photon character. Gen- 
eral properties of the emission from collective atomic 
many body states and their mapping onto the radiation 
field has been considered in [10| , while the mode structure 
for thermal ensembles of atoms was discussed in [Til ]. 

The role of radiative interactions between atoms has 
been studied since Dicke pointed out how the decay of 



collective excitations may lead to superradiant emission 
[l~2l . Other workers investigated temporal and spa- 
tial signatures of the emitted radiation pointing out that 
induced atom-atom coherence results in directional emis- 
sion of light, according to the shape of the atomic sam- 
ple [lj-[l8[. Both the limiting cases of single excitation 
[l5j and of a completely inverted medium [lj| have been 
treated. A review of the theory of collective spontaneous 
emission and related works can be found in [19j. 

The goal of controlling light sources at the single pho- 
ton level has recently put renewed focus on the treat- 
ment of the decay of a single collective excitation stored 
in an atomic gas [20h34I | . The collective Dicke state, a 
single excitation symmetrically shared among all of the 
atoms, decays with superradiant character only if it is 
stored in an atomic cloud with dimensions much smaller 
than the wavelength of the emitted radiation. The ex- 
citation remains trapped in the opposite limit of a large 
ensemble 30]. In this latter different quantum 

state, often called a symmetric timed (or phased) Dicke 
state [llj], shows fast decay in a given direction. In this 
work, we show that phase matched emission also occurs 
for multiply excited symmetric timed Dicke states in a 
large enough atomic ensemble. When two or more atoms 
are excited in Rydberg states, however, interactions de- 
phase the atomic synchronization imparted by the exter- 
nal laser fields. As a result their phase-matched emission 
is suppressed. The role of virtual processes, neglected 
in earlier works, has been shown to modify the decay 
of trapped states [3(| Hl| . Conditions and protocols re- 
quired to create such state have been recently investi- 
gated [2(], [24|. Furthermore, new systems have been ex- 
plored which show very interesting fight emission prop- 
erties in novel trapping geometries [HI, H^] . 

The paper is organized as follows. We briefly review 
the theory of interaction induced spin wave dephasing ac- 
cording to in Section |TT| this discussion motivates our 
focus on the phase-matching condition for multiple emis- 
sion processes and the influence of atomic interactions 
on it . In Section [ill! we discuss the Hamiltonian for the 
interactions between atom and field in the retrieval pro- 
cess. Section IIVI contains the main results of the treat- 
ment of the single excitation [22, H3 which is briefly re- 
viewed in Appendix |XJ The latter gives context for the 
analysis of multiple excitations in Section [V] (two exci- 
tations) and I VII (n excitations), leading to identification 
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of a phase matched radiation mode onto which the spin- 
wave properties are mapped. We summarize the results 
with final remarks in Section I VIII Appendix [B] provides 
technical details of the Wigner-Weisskopf approximation 
in the case of multiple atomic excitations. 



II. INTERACTION INDUCED DEPHASING OF 
MULTIPLE SPIN WAVES 

We consider a system of N atoms. We focus on three 
single atom energy levels that are coupled by light fields: 
a metastable (Rydberg) level \r), an intermediate level |e) 
and the ground level |<7), see Fig. [TJ The ground state of 
the atomic ensemble is the product state \G) = \gi....gpf); 
from here-on we will write atomic product states listing 
only those atoms excited out of the single-atom ground 
state, e.g., |s Ml ...s' ), where s...s' € {e,r} and the in- 
dices (/ii,. ../i n ) must all be different since the creation 
of two excitations on the same atom is forbidden. Tran- 
sitions between different atomic levels are described by 
the single-particle operators aff — |s /J )(s' |. We define 
collective excitations of level \r) in terms of spin waves, 
whose destruction operator is given by 
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(1) 



with r M the position of atom /i and ko the wavevector 
associated with the excitation. We also refer to this op- 
erator as the annihilation operator for symmetric atomic 
excitations that in Sections IIV1 fVl and I VII we relate with 
the symmetric states in the timed-Dicke basis. We con- 
sider atomic motion to be frozen, an approximation that 
requires the wavelength of the stored spin waves to be 
longer than the distance traveled by an atom during the 
storage and retrieval time Q. 

In order to describe the interaction induced dephas- 
ing of multiple spin waves, we consider a simple model 
with up to two excitations. Assume the initial excitation 
process brings the atomic ensemble to the state: 

|*o> = co\G) + Cl Sl\G} + C2 ^!|G). (2) 

Two-body interactions between atoms excited to level |r) , 
that take place after the excitation for a time T, lead to 
a phase shift 
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where = U^JT '/h is proportional to the two-body 
interaction strength . A signature of the multiparticle 
dephasing may be found in the two-particle spin wave 
correlation function defined by 
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(4) 



For a single excitation — 0. To illustrate the effect of 
dephasing, we calculate ((H) for the state resulting from 
([2]) after the phase shifts J3]) 



< 2 >(T) = 
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For @, a truncated coherent state c„ = 1/Ven!, 
n = 0, 1, 2, with $ M „ 0, 



7 ^)(0) = | < 1, 



(0) 



O...00, 



whereas for a coherent state, c n = 

g( 2 \0) = 1. For a random distribution of phase shifts, 
and a sufficiently long interval T > r, the sums in the 
numerator and denominator of ([5]) will vanish due to de- 
structive interference of the complex amplitudes. Since 
the denominator contains a constant term, we may take 
the approximation 
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The interactions can act to suppress the two-particle cor- 
relations as though only single excitations were present. 
The result is potentially useful as a source of single pho- 
tons if the spin-wave mode can be mapped onto a well- 
defined radiation field mode determined by ko, and if the 
dephased multiple excitations do not couple to this same 
mode in the radiative emission process. The following 
sections are devoted to analyzing these two issues, by 
identifying the states involved in the laser excitation, de- 
phasing and retrieval processes and their radiative decay 
channels. 



III. DESCRIPTION OF LIGHT-ATOMS 
COUPLING IN THE RETRIEVAL PROCESS 

The retrieval process is sketched in Fig. [T] A clas- 
sical laser pulse couples to the \r) ~ |e) transition and 
it is characterized by a Rabi frequency Ql, wavevector 
and angular frequency wi, while the states |e) and 
\g) are coupled to the continuum of quantized electro- 
magnetic (e.m.) modes. We label these modes with the 
index ip = (k, A), where k and A designate the wavevec- 
tor and the polarization, respectively; the energy of the 
photon is hujk = hck, where c is the speed of light 
and k = |k|. The general state of the field, containing 
to photons {4>\....<j) m ) in p different modes, is given by 

|*e.m.) = \<t>l-4m) = 4i 4ml°)/^ e Where 
the operator a\ is the creation operator for photons in 
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Figure 1. Schematic illustration of the retrieval of spin waves 
(wavevector ko) from a large cold atomic ensemble. The exci- 
tations are stored in the metastable state \r) (e.g. a Rydberg 
level): a laser with wavevector k^ is incident on the atomic 
cloud promoting atoms into the fast decaying state |e). The 
signal is detected according to a phase-matching condition in 
the direction given by ko — kt. 



mode 4> an d the factor e 0i---0"> = nf=i( n i') introduces 
the correct normalization depending on the populations 
of the different modes [Tij ]. The state |0) is the vacuum 
state defined by a^|0) = O,V0. 

The Hamiltonian for the coupled radiation field and 
atoms, in the eletric dipole approximation, is given by 
H = Ho + V, where: 

N 

H a = J2 E fiw«o£"+5^7iWka£ a *> (8) 

fi—l s—g,e,r 

ro N 

V = —Y J [e*^-'"-^*)*™ + h.C, 

AT 

"EE WkHH^ 1 " ) ■ (d eg n^;s) + h.c] . (9) 

We have defi ned the electric field per photon 8k = 
\J?ujJk/ {2e$V). The polarization vector for mode <j> is 
Erf, while d eg xi^ is the electric dipole matrix element for 



the transition |e)-|g) of atom ji. In the following, we in- 
troduce the radiation-matter coupling constant hg^^ — 
—£kd eg (e c f 1 ■ n M ). We make the rotating wave approx- 
imation (RWA) for both transitions. For the \r) — |e) 
coupling, RWA holds because we consider a resonant or 
quasi-resonant laser; for the quantized field, it has been 
shown that virtual processes have only a minor effect on 
fast deca yin g states, which are the focus of the present 
work [27l.l30|| . For this reason we make the RWA allowing 
to simplify the problem as the total number of excitations 
in the system, atomic plus photonic, is conserved. In the 
following sections, we consider separately the dynamics 
of states with a different number of excitations. 

IV. SINGLE EXCITATION 

The theory of a single excitation of the atomic ensem- 
ble is discussed in appendix [A] Here, we summa- 
rize the most relevant results. Using the notation pre- 
viously defined, the state vector may be written using a 
product state atoms- field basis as, 

|*> = E G |G; 4>) + 0) + 0)) • (10) 

4> M 

Applying the Wigner-Weisskopf (WW) approximation 
[37j | to the equations of motion (IA1|) . we derive a sys- 
tem of equations for the amplitudes of the excitations in 
level |e): 

^E = -M1 y" dt'e^n L {t')E,{t') 
F N 

-jEW- ( n ) 

!/=l 

where Y = uj 3 g d1 g / (Sneohc 3 ) is the single-particle decay 
rate for the excited level |e). We have defined Au>k — 
u! eg —LUk, Aul = cj re — ujl and r = t' — t. We define u) 8a > = 
ui s — uv . The function / is defined in (|A3[) according to 
ref. Jig. 

Employing the timed Dicke basis (|A4|) rather than the 
atomic product states as above, leads to equations for the 
amplitudes St defined in [A] In particular, for the sym- 
metric state (I = 0), the unique state in which all atoms 
share the excitation with equal probability, we obtain 
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^%Ko + E^ M1 (kK, 



(12) 



We have introduced coupling functions, V^'^k), i = matched direction k„ = k - k L , 
0,1,2,..., all of which depend implicitly on the phase 

e^k) = E^ (k ~ k ' oHr *~ M ' 



N 



V£*(k) = C NC £ e^'^ Y: Sf>(k - k' ) 



(13) 
(14) 
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Figure 2. (Colors online) Comparison of coupling strengths for the decay of a single spin wave into symmetric and non- 
symmetric modes from Eq. (|12[1 . Cubic sample, side L — 10/im, N = 100 particles (p = 10 11 cm -3 ). a) Coupling between a spin 
wave and one of the non-symmetric Dicke state (£ = 99): S = S^/ \/T{£+l) according to (|15|l . b) Normalized coupling 
of the spin wave to itself, Vgo'^k), see Eq. (|13|l . c) Comparison of normalized vj^^k) (colormap) with the total coupling 
to non-symmetric modes Y2e ^ol'^k) (red surface) : we assume equal amplitudes for all the modes. We assume k.' \\ z. The 
colormaps are based on the minimum (blue) and maximum (red) values of each plot. 



where 

S e p {k - k' ) = (V (k »- k ^ - e l ( k o- k )'^+i) . (15) 

The normalization coefficient is Cnc = 1/ y/Nl(t + 1). 
Some illustrations of the coupling functions are shown in 

Fig. m 

The expression (QUI) is peaked for k = k' and decreases 
far from this condition due to destructive interference of 
the different atomic phases. In the limit of a large number 
of atoms, it can be thought of as a non-zero width Dirac 
delta function. As pointed out in 27]. the coupling (IT41 is 
the product of two terms. The first sum is the same as in 
(fP5|) and is peaked in the vicinity of k = k^, while all the 
Sp functions are strongly suppressed in that region of the 
wavevector space, as it is evident from (1151) . The result 



is that in the equation of motion (fl"2)) the self-coupling of 
the symmetric excitation dominates over coupling to the 
other, non-symmetric, basis states. Therefore, we omit 
these terms and obtain an equation for the decay of the 
symmetric state with the generalized decay rate jl6j 

IV = r 1 1 + i E e-* k ° ( ">/> I . (16) 

\ M ^fJ, J 

This quantity is in general complex and contains both the 
superradiant broadening, Re(rjv), and (Lorentz-Lorenz) 
frequency shift of the atomic transition, Im(rjv). 

Recall that retrieval of the spin wave stored in the 
metastable (Rydberg) state |r) is done by means of a 
single-particle laser 7r-pulse, with wave vector k^ and du- 
ration T, followed by radiative decay from state |e). In 
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the case of a fast pulse, T <C l/Re(IV), the probability 
amplitude for field mode (f> is given by, 

p -(r N /2+iAu k )t _ i 
G0 W = ff^ ^ /0 I ■ a ^oc(k). (17) 



Tjv/2 + jAwfc 



We have defined 



(18) 



which is a Fourier component of the transition ampli- 
tude between the symmetric timed Dicke state |Eo) and 
the atomic ground state. The behavior of this term is 
illustrated in Fig. [3] Equations (fT7|) and ([T8|) define a 
phase matched mode [lol . [ill |22| into which the emission 
of light is concentrated. The non normalized creation 
operator for this phase matched mode is given by, 



Vo G (k) 



it 



(19) 



with a Lorentzian distribution of emission frequencies 
(fP7)) . of width r^r. The information stored in the spin 
wave is retrieved in the direction given by the phase- 
matching condition k = kp. We can summarize the exci- 
tation and retrieval in a simple sequence: 



rj«-,«(k-k' i 



Gr( k o.A)_ 



(20) 



DOUBLE EXCITATION 



Here, the normalization factors are chosen so that the 
summations run independently over all allowed values 
for the atomic and electromagnetic labels. We use the 
notation (y , fi) to indicate a sum on both indices with 
the exclusion of the term v = fi, in order to distinguish 
from the notation v ^ /i that indicates a sum over the 
first index, different from a given value of the second one. 
The equations of motion for the amplitudes are: 



dt 



8E* 
dt 



dt 



= J2 9<t> e ' lAuikt [e 4k ' rfl £* + e l ^Ef] (22a) 



(22b) 




-i\t-T u _ — iAujkt 



Ef' 



-g#e- ik "»e- i * w >>' t E* 



(22c) 



Since we assume that the system is initially loaded with 
two atomic excitations, we start our analysis from (|22a[) 
and (|22bl) . According to the discussion in appendix [Bj 
we may neglect the coupling terms depending on . 
We separate the process of emission of the photons in 
two steps: from the double atomic excitation to a single 
excitation and then to a two-photon state. The equations 
of motion for double atomic excitations are 



Before treating the general case of an arbitrary number 
n of excitations (where n << N), we study double exci- 
tations in detail. This contains the essential arguments 
of the general case shorn of the heavier book-keeping. 

The equations of motion for the three-level atomic 
structure with two excitations are much more compli- 
cated than the single excitation case. In principle, we 
should take into account all the possible combinations 
of two excitations. However, as discussed in the single 
excitation case, and in appendix [X[ it is possible to sep- 
arate the excitation of the metastable state |r) with the 
7r-pulse from the successive decay of the excitations from 
the state \e). Furthermore, the Hamiltonian for the laser 
coupling is symmetric with respect to different atoms and 
it preserves the symmetry of the initial state. With these 
observations, we may reduce the problem to the pair of 
states |e) and \g). 

We first write the doubly-excited state in the atomic 
product state basis as 



dE^ 
dt 



= "E4 fdt'e-^ 
Jo 



.J^e^-^E^t') 



J2e ik < r »- r ^E av (t') 
(23) 



This set of equations is the analogue for two excitations 
of (|A2[) for the single excitation: the amplitudes of all 
the possible pairs are coupled together through the con- 
tinuum of field modes. In this case, we have suppressed 
the coupling to the laser field, as discussed previously. 
A way to proceed at this point would be to apply the 
WW approach and then diagonalize the resulting system 
of equations. As we are interested in symmetric states, 
however, we choose to investigate the couplings between 
symmetric and non-symmetric collective atomic excita- 
tions, following the treatment of the single excitation. 



0,0' 



Km) 



0>m 



A. Timed Dicke basis for double excitation 

(21) 

We introduce the timed Dicke basis for the atomic 
states. For a single excitation we use (|A4I) . while for 
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the double excitations, the basis is given by: 



-oral. 



N 



e u e 



(24a) 



E ^l) = [ e4k " (rT(1,+ ^ <2)) | e 7(D e 7(2)> 



ik 0'( r *+l(l)+ r « + l(2)) 



(i)+r,+K2)i| ef+1(1)ef+1(2) ) . (24b) 



Here I G [1,W - 1], with TV = (^) . The labels f and 7 
are used to label pairs of atoms; l(V) and 1(2) indicates 
the first and second atoms in the pair; similarly for 7(1) 
and 7(2) [38[. The subscript [2] in square brackets refers 
instead to a two atom excitation state. 

We rewrite the state (f2"Tj) in the timed Dicke basis as 

I*) = £ -5- G**'\G;M/) + E E 4lE.;0)+ 



AT-1 



E <^[2]l E f[2];0). 



(25) 



£=0 



We derive the equations of motion for the new ampli- 
tudes. As above, we integrate the equation for Sf(t) and 
we substitute the result into the equation for <^ 2 ] after 

dropping the contribution from 



dS, 



[2] 



dt 



E^EE^"^ 1 ^!^ 

rf> £' ,] fJ,,v 

ft 

J -iAu t T. 



x (E>f \E e , [2] ) / dt'e-^ T £ m {t'). (26) 
Jo 

The structure of the coupling coefficients suggests sep- 
arated decay channels for symmetric and non-symmetric 
states. We define 

4 2 / 21 (k)=^e ik '^- r '')(E, [2] |^ar|E f[2] ), (27) 



and calculate the couplings involving at least one sym- 
metric state, 



^ [ o' 2l (k) = ^E el(k ^ Hr ^ r '^ (28) 
V™(k) ^K^e^^^tk-kJ), (29) 



7=1 



We have defined the normalization coefficient Cue 
1 / VAfC and the functions 



S^2i(k - k' ) = e-^-^-'TW + e -*(k-^)-'7(») 

2, ±1!^, 0. (30) 



> 7 [2]\ 

e -i(k-kJ)-r <+ i(i) 



-»(k-k{,)T/ +1(2 ) k - >k t), 



These expressions represent the leading order, while cor- 
rections are 0(1/N). This result shows that the symmet- 
ric excitations decouple from the non-symmetric states 
also in this case. In particular, Eq. (|2"5|) is sharply peaked 
at (5(k— kg), the detailed form depending on the atomic 
distribution. The functions (f30| vanish for k = k , caus- 
ing suppression of ([29)) close to that condition. The sepa- 
ration of symmetric and non-symmetric modes is related 
to a phase-matching condition as for the single excitation. 

We write ([26)1 for I = taking into account only the 
leading coupling, 



dt 



E~4 



^ e -,A^2,2] (kK ( ^ (31) 



From (|T6|) , we see that the amplitude for the double sym- 
metric excitation decays with radiative width 2Tn [22] . 
with the photon emitted in the phase matched direction 
determined by wavevector k . 

We consider the equation of motion for the amplitude 
of the single atomic excitation plus one photon: 



~dt 



e-^{^\£ m )S m . (32) 



The coupling between the singly and doubly excited 
atomic states is fixed by the couplings 



4 1 ; 21 (k) = ^e- <k ^(E,|af|^ 



[2]/- 



(33) 



Among these matrix elements, we calculate those that 
contain at least one symmetric state, 



V& 2] (k)=C NC J2 S ~ l l2 ] (k-K), (35) 

7=1 

1 

v£ ,21 (k) = Cmc E i-siO* <)] • ( 36 ) 

These expressions determine how a double atomic exci- 
tation decays to a single excitation with emission of a 
photon. The dominant term is (|34|). which is peaked in 
the preferred direction k = k . The phase matched direc- 
tion favors the coupling between the symmetric excita- 
tions with respectively two and one excited atoms. This 
decay channel is separated from non-symmetric atomic 
states: the S functions (fT5[) and (I3TJ1) . for single and dou- 
ble excitations, appear in (|36|) and (|35|). suppressing the 



contribution of these states to the emission around k . 
We specialize (|32"T) to the case £' = corresponding to 
a double spin wave and reintroduce the coupling to the 
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two-photon continua, giving 



dS[ 
dt 



9<pe 



-•AUfctT/I 1 . 2 ] 



var J (kK0[3]- 



(37) 



The matrix elements Via couples single atomic excita- 
tions and atomic ground state; we have already calcu- 
lated the form for the symmetric state Vqg in (|TH|> . The 
transition amplitudes to non-symmetric states are, 



F, G (k') = -^X>i(k'-k' ). 

V a— 1 



(38) 



These results confirm the behavior found for a single 
excitation and show how the radiative interactions among 
atoms drive the decay of two spin waves (doubly excited 
symmetric timed Dicke state). 

We can summarize the emission processes by the se- 
quence: 

v o r~'(k-ko), ^,A) y, ~a(k'-kj) ; GlK , x)lM 

(39) 

All other couplings are proportional to some function 
/(k - k' ) -> for k -> k[,. 



discarding interaction terms between , 
cussed in appendix |Bj We obtain 



and 



as dis- 



dSf 
dt 



-^ e - iA "<=Vi 1,2I (k)# 



[2] 



- E 4'^c(k')^(k') T dt'e- iA »^4{t'). (43) 

[1 21 

Recall that the coupling V2, ' is suppressed with respect 



in, and T^g is dominated by Vqg- We solve for to 





the leading order. The WW approximation gives 



(-iAco k -T N )t 



■v r /9 — ^b[2](0). (44) 
lAujk + 1 at/2 1 ' 

In the large ensemble limit the emission is strongly 
peaked in the phase-matched direction as previously 
stressed. 

Finally, we solve for the two-photon amplitude . 
In the asymptotic limit t — > oo, we obtain 



1 



5^o [ o' 2l (k) 



iAuk + Tat/2 iAujk' + Tat/2 



^0[2](0). 

(45) 

By comparison with (ITT)) , it is clear that the decay of the 
double spin wave results in a pair of photons emitted into 
the phase matched mode, (fT9")l : 



I*) 



^(blf\G;0). 



(46) 



B. Decay of a double spin wave 

Based on the approximations exploited above, we give 
the explicit temporal evolution of the amplitudes in- 
volved in the decay of a double spin wave. The initial 
conditions are: 



(40) 



By switching to the timed Dicke basis and through ap- 
plication of the WW approximation, the solution of (|3"Tj) 
is 



^0[2](*) 



-r N t 



<£b[2](0), 



(41) 



where Tat is defined in (1161) . In order to solve for and 
G^', we use (|37J and 



" f 'dt'y2\g 4 e- iA ^'v Ge (k)4'(t') 
Jo t l 

+^ e - iA ^'*V G ,(k')<(t')}- (42) 



We substitute this expression into (|37|) . retaining the cou- 
pling between states with different single atomic exci- 
tations and the same photonic mode, St and St , but 



Effect of atomic interactions 



We finally consider the effect of the phase shifts 
§S§ induced by two-body interactions for atoms in the 
metastable state \r). The state | E [2] ) evolves into 

I*) = / E e ik °'^ +r ^*-M. (47) 

v—l ii>f 

We expand |$) in terms of the symmetric and non- 
symmetric modes: 

(WH^E^ 



-t[2] 



**j(lb(2) 



£+1(1)4 + 1(2) 



(48) 



(49) 



3=1 



In the limit — > 0, we obtain the previous results: the 
amplitude for the symmetric timed Dicke state is 1 while 
all the other terms vanish. In the case that the phase 
shifts are large and broadly distributed, the amplitude of 
the symmetric state is quenched indicating the emission 
into the phase matched mode. If we recall the state ©, 
this quenching of the two-photon amplitude in the phase 
matched mode corresponds to the suppression of the spin 
wave correlation function, consistent with single photon 
emission in direction kn. 
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Figure 3. (Colors online) Behavior of the normalized emission amplitude VoG(k) in the phase matched mode. We consider 
k = (k,0, \Jk' 2 — k 2 ) and kp = (0,0, k' ). a) Effect of diffraction from spherical ensembles with different radius 7? at fixed 
density p — 10 12 cm -3 . Green dot-dashed line R — 20pm, Red solid line R = 10pm, Blue dashed line R — 5pm. b) Dependence 
of the peak on the density of the atoms, for a spherical cloud, R — 10pm. Red solid line p = 10 12 cm -3 , Green dashed line 
p — f0 11 cm -3 , Blue dot-dashed line p = 10 10 cm -3 ; for the former two cases the peak is completely symmetric with respect x 
and y directions while for the latter density, we also plot the case k = (0, k, \Jk2 — k 2 ) which is the Black dotted line. This 
shows the effect of fluctuations at low density, c) Signature of the shape of the atomic cloud in the tails of the phase matched 
emission. Density is fixed p = 10 12 cm -3 . Red solid line is a spherical ensemble, R — 10pm, Blue dashed line is a cubic 
ensemble, side L — 20pm, and Green dot-dashed line is a spherically symmetric Gaussian distribution with <J x , y ,z = 2.5/xm. 



VI. MULTIPLY EXCITED STATE 

In this Section, we show that the results of Sec. IIVI 
and [V] hold in the more general case of n <C N spin 
waves stored in the ensemble. As for the double exci- 
tations, we assume that a laser pulse transfers the spin 
wave excitations from |r) to |e) before the emission takes 
place. In order to proceed, we refine the notation: we 
label the many-body atomic state through the indices of 



excited atoms |e, 



\pLl...pL n ). 



We consider the symmetric initial condition 



^i...p»(0) 



Jk' -(r„ 1 +...+r^ Il ) 



(50) 



while all the other amplitudes vanish at t = 0. 

We introduce the timed Dicke basis for m = l,....,n 
material excitations: 



-Ofrol 



1 J_ V e< 

777 ' * ' 



(r M1 +...+r M 



(51a) 

(r7«+-+ r T(m))| 7 (l)... 7 ( m )) 



_ e a&-('ww+™+»«+Km))|^ + i(i)..^ + i( m )) . 

(51b) 

Here, following the case of double excitations, we define 
£ and 7 as indices for collections of m different atoms, 



1 e 

E W=^E[ elk ° 



£(l),£(2),...,£(m). The notation [m] is a reminder that 
the state contains m atomic excitations. 

By using this basis and similar approximations to those 
described for the treatment of double excitations, we ob- 
tain an equation for the decay of the nth-excitation state, 



(:{„] 



dt 



E 4 E V U>' n] W /* dt'e-^S t[n] (t>), 

(52) 



which has the same structure as for the double exci- 
tation. Here, we define the radiative coupling among the 
timed Dicke states by, 

V#> n] (k) = e ik ' (r "- r " ) (E, W |^af \E e , [n] ), (53) 

giving explictly, 

^(k) J2 e^H'^ + O fl) , (54) 



4"'" ] (k) = ^=E{E el(k - H) ' r "^N( k - k o) 



(55) 



where we have generalized the definition of the functions 
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(HBl and (Ei) to 



^[„](k-k' o ; 



O *(k-ko)-r 3(s) _ i(k-k^)-r f+1(a) 



s=l 



(56) 

These expressions make evident also in this case the sep- 
aration of the decay channels between the collective sym- 
metric excitations and the rest of non-symmetric atomic 
states, as shown in Fig. [2] for the case n = 1. In fact, 
Eq. (fMf has a sharp peak for k = kg where the couplings 
in (|55|) all vanish. We obtain from (|52|) that the symmet- 
ric timed Dicke state | E [n] ) decays with rate nRe^Ar). 

By following the calculation for the double excitation 
case, we evaluate the decay process which leads from n 
excited atoms to n — 1 plus one photon, in order to char- 
acterize the emitted radiation. By analogy with Eq. (|32|) , 
we define the matrix elements connecting the symmetric 
state | E "' ) to the basis with n — 1 atomic excitations, 



Those couplings involving at least one symmetric state 
are, 



(58) 



vi B -^(k) = J £ S* [n] (k k), (59) 



/iVT7E[-4n-l]( k 0- k )]- ( 6 °) 
\n) L 7=1 



These have similar properties to (|34[) . (|35|) and (|36l) . In 
fact, (155)) is peaked for k = k showing the predominant 
emission in the phase matched direction where the cou- 
plings (|59"|) . ([TO]) are strongly suppressed because of the 
3\ n \ mnc ^ ons - This indicates the first photon in the de- 
cay of the multiple spin wave will be emitted in direction 
k and the atomic state will retain its collective symme- 
try. This result for the n —> (n—1) de-excitation, together 
with the results for the single and double excitations af- 
firms by induction that the subspaces of symmetric and 
non-symmetric atomic states are separated in the radia- 
tive coupling with the quantized electromagnetic field. 
The argument is valid for any number of atomic exci- 
tations throughout the whole decay process. The final 
state of the system is then given by 



I*) 



(&L) n |G;0). 



(61) 



We stress that the calculation is valid only for a large 
number of atoms since we have neglected corrections of 
0(n/N). 

Two-body interactions that dephase the symmetric n 
atom excitations in the metastable state will once again 



act to quench emission into the phase matched mode. 
The state resulting from the dephasing of n spin waves, 
with pairwise phase shift &fj, u , is given by 



V It!/ fJ,\> . . .> fJ. n 

(62) 

The amplitude of the phase matched symmetric state is 



(E 0[n] |$ (w) ) 



N 



e 1 ^*^, (63) 



fj,i>...>p,„ 



which vanishes in the limit of large and broadly dis- 
tributed phase shifts, except for the single excitation 
amplitude that results in single photon emission in the 
phase-matched direction. 



VII. CONCLUSION 

We have described the radiative retrieval process of 
multiple spin waves stored in metastable Rydberg states 
of an atomic ensemble. We have shown that the de- 
cay from a weakly excited ensemble is strongly di- 
rectional thanks to the enhanced coupling to a phase 
matched mode and suppression of the contribution of 
non-symmetric excitations. The demonstration given is 
valid for an arbitrary number of excitations much smaller 
than the total number of atoms. The use of the timed 
Dicke basis allows to focus on the relevant phase matched 
mode and can be expressed in terms of the transition am- 
plitudes between the timed Dicke states mediated by the 
electromagnetic field. As expected, the radiative cou- 
pling to the phase matched mode contains information 
about the size of the ensemble and its fluctuations are re- 
lated to the density and the shape of the cloud. The anal- 
ysis enables a mapping of atomic excitations into emitted 
photons. It indicates that symmetric atomic excitations, 
those created by laser driving, undergo phase matched 
emission. A pairwise interaction of metastable stor- 
age atoms, however, suppresses these symmetric ampli- 
tudes and quenches multiphoton emission into the phase 
matched mode, leaving only the unperturbed single pho- 
ton emission process. This insight may have application 
to fast single photon sources based on cold atomic en- 
sembles. 

We acknowledge financial support from NSF and 
AFOSR. We thank Y. Dudin, A. Kuzmich, H.H. Jen and 
S.D. Jenkins for useful discussions. 



Appendix A: Review of the theory for a single 
excitation 

The case of a single excitation stored into an atomic gas 
has received extensive attention in the recent years. Here, 
we limit the discussion to review the relevant results 12 21, 
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|27I | in a way appropriate for generalization to multiple 
excitations. We also use this example to show how to 
eliminate the dynamics of the metastable state \r), thus 
reducing significantly the complexity of the problem for 
the case of multiple excitations. If a single excitation is 
present in the system, the general state may be written 
as (fTQ|) . The equations of motion for the amplitudes of 
the different basis states are 

dC^ 

ih— = -mJ^a^e-^e-^E^, (Ala) 



ih 



dt 

dE^ _ hfl L _ ik 



g-iKL-rpg-iAwitn 

2 M 

-ih^g^-^e^'G^, (Alb) 



dt 



dt 2 M 



(Ale) 



Here, we introduce the notation Awj. = u> eg — uj^ and 
Aojl = oJre — Wi. We integrate (|Alc[) and (jAlal) and we 
substitute the results in (IAlb[) : 



dE M = n L {t) 

dt 4 



dt'e^^Lit^E^t') 



- 4 f dt'e- iA "» r £ e ^-^E v {t'). (A2) 

A. JO 



1. Timed Dicke basis for single excitation 



The timed Dicke basis for a single excitation in the 
ensemble is defined as follows I22| 



1 N 



(A4a) 



l R ^> = 77 E [ elk0 ' rM M - e ^°' re+1 l s ^+i)] i ( A4b ) 

V u=l 



herefe [1,JV-1], £ = £(£+1). If the excitation is in the 
level |e), we substitute R —¥ E and we use the wavevector 
k' = k — ki . In general we refer to the I = state as the 
symmetric state since the excitation is shared between all 
the atoms with equal probability, while all the other basis 
vectors are non-symmetric in this sense. We rewrite the 
state (TTU)) in terms of this new basis 

N-x 

= ]Tg*|G;^) + («|E 4 ;0)+« / |R<;0)). (A5) 

A 1=0 



We write the equations for the amplitudes of the new 
basis states: 



Here, we assume that the coupling coefficient is 
independent from [i, which is equivalent to consider 
all the atomic dipoles aligned in a given direction; 
this model can be realized by properly choosing the 
external fields used for excitation and trapping of the 
atoms in order to allow for a single transition to be 
active. We also define r = t' — t. We substitute the 
sum over the e.m. modes with an integral, Yl<f>=(k \) ~~ ^ 
Ea=i,2 V/(2ir 3 ) da /* d6 sin 6 <W/c 3 , where 
the angles (a, 9) characterize the direction of the 
wavevector of the photon and A labels its polarization; 
the sum over the polarization for the coupling constant 
gives Ea=1,2 n 9% ^i 1 ~ cos2 °) dlJ( 2e oV). By using 
the Wigner-Weisskopf approach [33], we obtain (fTT]) 
from (|A2I) . In this latter equation, the function / is 
defined as [l6( 



271 



f* v = —j da I d9sin6{l-cos 2 6)e 



(A3) 

and it depends on the relative distance between particles 
fi and v as well as on the wavevector of the emitted pho- 
ton. Equation (fTTI) corresponds to the expression derived 
in literature |l6l . l30ll3"lj ]: it is possible to analytically solve 
it in terms of Bessel functions for the case of a scalar pho- 
ton. Since we are interested in the decay of a spin wave, 
we proceed in a different way by introducing the timed 
Dicke basis for the atomic wavefunction. 



dG^ 
~dt 

dt 



d^ e 
~df 



9*£ 



ikr. 



{G\df\Et)g t , (A6a) 



Y g4,e lAuJkt e lk -^(E e \af l 9 \G)G' l> , (A6b) 

H,6 



2 1 



(A6c) 



We integrate (|A6a[) and substitute the formal result for 
G^ into (|A6bl) to determine the coupling between the 
different timed Dicke states mediated by the field. We 
analyze the behavior of the transition amplitudes for en- 
semble emission and re-absorption between the symmet- 
ric timed Dicke state £ = with itself, £' = 0, and with 
the rest of the basis, £' ^ 0: 



(E |^|G)(G|<7f |E ) 



,iko-(r„-r M ) 



N 



(A7) 



1 N 

(E \a^\G)(G\a^\E e7t0 ) = -y=J2 



-iko-r Q r 



e lk ^ +1 W . (A8) 



By using these matrix elements, we are lead to (fi 
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2. Retrieval of a single spin wave 



As initial condition, we assume a timed Dicke state 
stored in the metastable state with wavevector ko 



i? M (0) = ^=e lkoI \ (A9) 



while all the other amplitudes vanish at the intial time. 
We may model the retrieval of this spin wave by means 
of a 7r-pulse from the laser field f2 l followed by the decay 
from the state |e). We consider a resonant pulse, lul — 
io re . Furthermore, we assume a pulse whose duration 
T ss 7r/£li, where Ql is the temporal average of the 
laser Rabi frequency, is much smaller than the ensemble 
decay rate, T <C l/Rc(TAr). Within these conditions and 
neglecting the weak coupling to non-symmetric states, 
the solution of (fT2j) may be written as 



calculate the derivative in time of the function (|A10|) : 

d£ n L (t) r* . ^ .^ 
-dt = -—j dt —^) 

(All) 

where we have used the fact that d t f3 — flL{t)/2. We 
remark that the approximation holds for 

^ T ~f[ dt'A cos W)\ e'^' « T ~fT « 1, 

(A12) 

which is exactly the regime we are assuming. It is also 
possible to allow for a temporal dependence of the coef- 
ficient A which provides an integro-differential equation 
for this quantity which turns out to be constant in the 
case (|A12[) . The temporal evolution is then separated in 
two parts: for t % T, the amplitude for the symmetric 
excitation in the intermediate state £q grows following 
the laser pulse, while for t ^> T the coupling to the con- 
tinuum of the radiation modes transfers the excitation 
to the field. In this temporal range, we can approximate 
sin(/3) 1 and the wavefunction asymptotically coincides 
with the e.m. field amplitude (fTTj) . 



£'o(t)=A 



n L (t') 



dt' 



(AlO) 



the factor A is set by the intial condition A = ^o(O) = 
1. We define the phase f3{t) = L dt'£lL{t')/2, and we 



Appendix B: Wigner-Weisskopf approach for 
multiple excitations 

We discuss the two approximations used in the applica- 
tion of WW theory to multiple excitations. We consider 
the case of double excitations, since it contains all the 
ingredients necessary to this demonstration. In order to 
show the validity of the approximations, we use the set 
of equations (|22f . We first focus on the coupling between 
the continua of one- and two-photon states. We integrate 
Eq. (|22cp and we substitute the result in (|22b|) to obtain: 



dE* 
dt 



T E lllJ - 



+94>'9<P e 



f t N 

, / ^^(t-Ore-a 1 '!',.-'.)^ 



(Bl) 



In the spirit of WW, we substitute the sum over the 
e.m. modes (j)' with an integral. This integral can be 
separated in a frequency and an angular part: the former 
one is used to perform the Markov approximation. We 
see the difference between the two terms in the bracket. 
In fact, the coupling to Ef depends on the phase factor 



e iAw fc ,(t t ) anc ^ together with the integral over ui' = ck' , 
it gives a function 6(t — t'). The coupling to the ampli- 
tudes E^ contains instead e lAWfc/t that results in S(t): in 
the solution of the equation, the contribution of this term 
vanishes. We are thus allowed to discard the couplings 
E$ -R> Ef . After application of WW, we are left with 
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the equation 



dEt 
dt 



■ e -i^{t-t') Ellv (t) 



V N 

-Y 

2 ^ 



(B2) 



where T is the single atom decay rate and we use the 
definition of the / function (|A3I) . 

A similar argument is valid to check the second approx- 
imation which consists of neglecting the terms depending 

I 



I 

on G**' when we substitute (I22b|) in (j2"2"aj) . We consider 
the expression in (|B2|) . We formally integrate the expres- 
sion and we substitute the result in (|22al) . which yields 



dE, 



)J,V 



dt 



£4 fdt> 

j. Jo 



Ys^ kt 



e -iAw fe (t-t') e -i(i-t') 



dt ' e -i(t-*')E 



(B3) 



We notice that the term in the second line of this equa- 
tion contains the same type of sum over the continuum 
of radiation modes as the last term in (|B1[) . In fact, the 
integral over the modulus k together with the oscillating 
term e lAwfc (* _t > gives a 5(t — t') in the terms in the first 
line, while its combination with e lA " fct in the second line 
produces the function S(t): these latter terms then van- 
ish when we perform the Markov approximation within 
WW and they do not affect the decay of the double ex- 
citations. It is possible to understand this effect in the 
following picture: the linewidth acquired by the states 
E^ as a consequence of the coupling to the amplitudes 



is much smaller than the bandwidth of their contin- 
uum of states, responsible for the decay of E^ v . While 
studying the behavior of E^, it is then possible to dis- 
card the terms depending on G^ : this leads to Eq. ([23|) 
in the main text. Furthermore, it is straightforward to 
verify that these approximations allow to find the correct 
results for the case of N — 2 atoms and n = 2 excitations, 
which agree with the results shown in [l4[. These argu- 
ments extend in the same way to the multiply excited 
states. 
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